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Fermions in an optical lattice near a wide Feshbach resonance are expected to be described by an 
effective Hamiltonian of the general Hubbard model with particle-assisted tunneling rates resulting 
from the strong atomic interaction [Phys. Rev. Lett. 95, 243202 (2005)]. Here, we propose a 
scheme to unambiguously test the predictions of this effective Hamiltonian through manipulation of 
ultracold atoms in an inhomogeneous optical superlattice. The structure of the low-energy Hilbert 
space as well as the particle assisted tunneling rates can be inferred from measurements of the 
time-of-flight images. 



Fermions in an optical lattice near a wide Feshbach 
resonance provide one of the most complicated systems 
for ultracold atoms 0, 0, i, i, S, i, 0, S, Q . The strong 
atomic interaction induced by Feshbach resonance pop- 
ulates many excited lattice bands and causes direct cou- 
plings between neighboring sites. To understand this im- 
portant system, one needs to have a theoretical model. 
An effective Hamiltonian has been proposed in Refs. 

based on the arguments of the low-energy Hilbert 
space structure, which offers a significant simplification 
for description of this system. The effective Hamiltonian 
takes the form of the general Hubbard model (GHM), 
where the effects of the multiband populations and the 
direct neighboring interactions are incorporated through 
the particle-assisted tunneling rates, ft is important to 
test this model by comparing its predictions with experi- 
mental observations. However, such a comparison is usu- 
ally difficult because of the lack of exact solutions to the 
GHM and complications in real experimental configura- 
tions (such as the inhomogeneity due to the global trap) . 

In this paper, we propose an experimental scheme to 
quantitatively test the predictions of this effective model 
by manipulating strongly interacting atoms in an optical 
superlattice. The optical superlattice provides a pow- 
erful tool, which has been used in recent experiments 
for demonstration of the spin super-exchange interaction 
[Tol [Til . Il2| . In the experimental configuration with an 
inhomogeneous optical superlattice, through manipula- 
tion of the lattice barrier and the external magnetic field, 
we show that one can reconstruct the two-site dynamics 
from the measured time-of-flight images. The measured 
dynamics can then be compared with the exact prediction 
from the general Hubbard model, offering an unambigu- 
ous testbed for this complicated system. The proposed 
measurement also allows a complete empirical determi- 
nation of all the parameters in the effective GHM. 

We consider two component fermions (denoted by spin 
a = T: I) m an optical lattice near a wide Feshbach reso- 
nance. Although in general many lattice bands get pop- 
ulated due to the strong atomic interaction, we note in 
Refs. [5, 9] that for this system the low energy states at 
each site are still restricted to only four possibilities: ei- 
ther a vacuum denoted by |0), or a single atom with spin- 
<7 denoted by a J. |0), or a dressed molecule in the ground 



state \d) which consists of superpositions of two-atom 
states distributed over many bands. All the other states 
(such as the three-atom states or the dressed molecule 
excited states) are well separated in energy, and there- 
fore not relevant for low-temperature physics. Based on 
this low-energy Hilbert space structure and general sym- 
metry arguments, it is shown in Q that the effective 
Hamiltonian takes the form of the GHM: 

H = K u / 2 ) n i fa ~ 1) - Mi^i] (1) 

i 

where Ui = ^ a a\ a a ia , = a\-a m {a =|, \ for a =|, j), 
U characterizes the effective on-site interaction (defined 
as the energy shift of \d) with respect to the two-atom 
state on different sites), /Uj is the chemical potential (we 
keep its dependence on the site i for convenience of the 
following discussion, where a global trap induces a site 
dependent energy shift), t a is the conventional single- 
atom tunneling rate, and g\ and 32 denote the addi- 
tional tunneling assisted with spin-a atoms (those two 
terms come from the multi-band populations in the \d) 
state and the direct neighboring atomic interaction in the 
lattice [5j). In this derivation, we have mathematically 
mapped \d) to the double occupation state a^a^ |0) 1 
(though their physical compositions are different). 

To sec whether the effective GHM gives a good de- 
scription of the low-temperature physics for this sys- 
tem, it is important to test the predictions of the GHM 
through experiments. To have an unambiguous test, 
it is better to design a configuration such that the 
GHM allows exact solutions. The optical superlattice 
provides such an opportunity. To produce the opti- 
cal superlattice, one simply adds a 3-dimensional lat- 
tice V 2 = V20 H a =x,y,z sm 2 (7ra/2a) to a lattice V\ = 
V\q sin 2 (ttz I a — ip) in one spatial direction (say z) , where 
the periodicity 2a of V 2 is twice that of V\ pjl EH, E3] • 
If V\q is sufficiently large relative to V20, the superposi- 
tion of these two potentials produces a series of double 
wells along the z direction. The dynamics in each dou- 
ble well are independent of the others provided the bar- 
rier between wells (controlled by V20) is sufficiently large. 
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FIG. 1: The time sequences for the magnetic field (B) and the 
lattice potentials (Vio and V20) to achieve state preparation, 
controlled dynamics, and detection. (Note that the optical 
barriers are high during sweeping of the B field.) 

Taking the relative phase ip to be nonzero introduces an 
energy bias 612 between the minima of each double well. 

It is easy to calculate the dynamics in each double well 
from the GHM. However, it is unclear how to directly 
measure the dynamics without individual addressing of 
each well inside the lattice. The conventional time-of- 
fllght (ToF) images involve averages over all the potential 
wells. These signals are further complicated by the pres- 
ence of a global harmonic trap V g = J2 a =x y % ww^a 2 /2 
inevitable in an optical lattice, which makes each dou- 
ble well slightly different. In the following, we show a 
scheme that can map out the detailed dynamics in each 
double well from the measured ToF images even with the 
presence of these complications. 

The scheme here combines the control of both the op- 
tical potentials and the magnetic field (see Fig. [1] for 
illustration). First, we need to load each double well 
with a filling pattern that sets the initial condition of the 
dynamics. This is achieved at the BCS limit of the reso- 
nance. In this limit, the atoms are free fermions, and we 
can control the filling pattern by choosing the total num- 
ber N = A T +N l and the polarization P = (N^-N^/N. 
Then, we turn off all the inter-well dynamics by rais- 
ing the barrier (controlled by Vio and V20) and sweep 
the magnetic field to the unitarity region. The sweep- 
ing speed v is fast compared with the inter-well coupling 
rate but small compared with the lattice gap of V\ so 
that the levels in each single well adiabatically evolve. 
Near unitarity, we turn on the inter-well dynamics for a 
duration t by adjusting Vio to lower the central barrier 
of each double well. These dynamics give information 
on the underlying strongly interacting Hamiltonian. To 
determine the final state after the dynamics, the central 
barrier is raised again, and the magnetic field is swept 
to the BEC limit with a speed similar to v. Depend- 
ing on the particle number in each well, we have atoms 
or molecules or their mixture with negligible interaction 
at the BEC limit. The ToF images for those atoms or 
molecules are then detected to determine the final state 
after the dynamics during time t. 

To test the prediction of particle-assisted tunneling, 
we need to compare the free-atom hopping rate t a with 



t a 2 =t a + gi and t a3 =t a + 2gx + g 2 , where t a2 and t a3 
correspond respectively to the hopping rates of a spin-f 
atom from the site i to j when there is a spin- J, atom 
on one site or on both sites. Let us first look at how to 
measure the free-atom hopping rate t a in the Hamiltonian 
|T]) . For that purpose, we need one atom per double well. 
By choosing the polarization P = 1 and Vio = (so we 
have at this stage single wells rather than double wells), 
the equilibrium distribution of the free fermions at the 
BCS limit automatically gives this configuration. The 
total atom number N within the global harmonic trap 
V g needs to be below 

N max = (4tt/3) {E bg /2muo V) 3/2 , (2) 

where E^ g = 2y / Via7r 2 h 2 /8ma 2 is the band gap for the 
lattice V2 and we have assumed oj x — oj v — to z = lu. 
Then, one can adiabatically raise the potential Vio with 
a bias S12 so that the atom sits on the left-side well in 
each double well. After raising Vio, #12 is reduced to 
zero. The system is then moved to the resonance region, 
and after turn-on of the dynamics for a duration t, the 
difference between the fraction of atoms in the left-side 
and the right-side wells over the whole harmonic trap is 
given by 

Nl-N r 2 W *° V • Vo a (*\ 

where flu = ^Af + 4^/2^, A, ~ mw 2 azi. The summa- 
tion of i in Eq. ^ is over all the occupied double wells 
in the global harmonic trap (with Zi the z-coordinatc 
of the center of the double well), and each double well 
has a slightly different bias A, due to the trap potential 
V g . After the dynamics, in order to measure the popula- 
tions Nl and Nr, the atom of the right-side well can be 
dumped into an excited vibrational state (corresponding 
to the second band) of the left-side well by rapidly raising 
the potential minimum of the right well relative to the 
left (through control of the phase ip) 0, EL 03 ■ The 
populations in different bands are then mapped out in 
the BEC limit through measurement of the momentum 
distribution of free atoms with ToF imaging. From the 
measured populations Nl and Nr, one can easily deter- 
mine the tunneling rate t a . Fig. [U (a) shows the typical 
time evolution of Nl — Nr from the dynamics, for which 
the oscillation period is determined by t a and the damp- 
ing is due to the inhomogeneity from the global trap [13| . 
In the frequency domain, the signal peaks at 2t a , and the 
inhomogeneity causes many smaller peaks at frequencies 
above that of the dominant peak (see Fig. [2](b)). 

To measure the particle-assisted tunneling rate t a 2, we 
need two atoms per double well, one spin-| and one spin- 
|. This can be achieved with the equilibrium distribu- 
tion of free fermions at the BCS limit by choosing PwO, 
Vio = 0, and the total atom number N < 2N max (with 
N m ax defined in Eq. [5]). The double well is still turned 
on with a bias so that both of the atoms are prepared in 
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FIG. 2: Population difference between the left and right wells 
(Nl — Nr)/N for the case of one atom per double well, (a): 
Population difference vs. time (in the unit of h/t a ). (b): 
Fourier transform of the population difference (frequency in 
the unit of t a /h), calculated for a time duration of 20h/t a 
to give a frequency resolution of t a /20h. The peak occurs 
at a frequency v — 2t a /h. Inset: With the time duration 
increased to give a frequency resolution of t a /2000h, we see 
that there are actually many peaks, corresponding to the dif- 
ferent frequencies Qu- fin depends on the z-coordinate and 
thus each peak corresponds to a different slice of double wells 
parallel to the z-axis. The slices containing the most occu- 
pied double wells are closest to z = 0. That is why those 
peaks (which have the smallest £lu) dominate. Because t a 
can be determined from the dominant peak, it is not nec- 
essary to resolve the other smaller peaks. In calculation of 
the inhomogeneity effect, we assume a spherical distribution 
with a diameter of 30 occupied double wells, and take the 
following typical values for the parameters: t a — h x 170 Hz, 
m = 6.64 x 10 - 26 kg (for 40 K), ui = 2tv x 80 Hz, and 2a = 765 
nm. 



the left-side well. For the dynamics near resonance with 
the Hamiltonian |1| , the state at any time involves a su- 
perposition of three components: a double occupation of 
the left or the right well, and a singlet state of two atoms 
over the two wells. We can determine t a2 as well as the 
on-site interaction energy U from the difference between 
the overall fractions of double occupation of the left wells 
and of the right wells, (N 2L - N 2R )/(N 2L + N 2 r). (Here 
N 2 l and N 2R are the total number of double wells in 
which the left and right wells, respectively, are doubly 
occupied.) These fractions can be directly measured at 
the BEC limit, where the double occupation of a site is 
mapped to a molecule state, and the molecules in the left 
and the right wells are distinguished through the band 
mapping and the measurement of the momentum dis- 
tribution (similar to the discussion above for the atomic 
case) . The single-atom occupation of a well is mapped to 
an atomic state at the BEC limit. Because of the large 
detuning between the atomic and the molecular state, 
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FIG. 3: The difference between the fractions of doubly occu- 
pied left and right wells (N 2L - N 2 r)/(N 2 l + N 2 r) for the 
case of two atoms per double well, (a): Population differ- 
ence vs. time (in the unit of h/t a2 ). (b): Fourier transform 
of the population difference (the frequency resolution is 1/20 
in the unit of t a2 /h. The peaks occur at a frequencies v\ — 

(yjP + lQt 2 a2 - U\ /2h and v 2 = (^/U 2 + 16t 2 a2 + £/) /2ft 

(in the figure we take U = 3t a2 as an example). Increas- 
ing the frequency resolution would reveal a series of smaller 
peaks on the high frequency side of the large peaks (as in Fig. 
[2]), but it is not necessary to resolve these smaller peaks to 
determine t a2 and U. 

the atomic population does not contribute to the time- 
of-flight imaging signal of the molecular fraction. 

The typical time evolution of (N 2L -N 2R )/(N 2L +N 2R ) 
is shown in Fig. [3] (a). In the frequency domain (see Fig. 
[3] (b)), one can see two distinct primary peaks in the 

Fourier transform, centered at (^\/U 2 + I6t 2 2 ± UJ /2. 

The smaller peaks from the inhomogeneity of the global 
harmonic trap do not obscure these two dominant peaks. 
The frequencies at which these two peaks occur can be 
understood by the fact that while the oscillation fre- 
quency varies from well to well due to the z-dependent 
bias, there are a greater number of occupied double wells 
near z = than for any other z-coordinate. Thus, the 
dominant peaks correspond to the zero bias case, where 
we have 

% 2 * = ^ cos (n_t) + ^ cos (n+t) , (4) 

with Q.± = (m 2 i ± U) /2ft and il 2l = ^U 2 + 16t 2 a2 /h. 

One can also measure the parameter g 2 in the Hamil- 
tonian (Jl]), which requires three atoms per double wells 
(two spin-|, one spin-|). One can consider this case as a 
single spin-f hole in each double wells, with a hole hop- 
ping rate of t a 3- This hopping rate can be measured by 
the same method as for measurement of the free atom 
hopping rate t a \. To prepare three atoms per double 
wells, one can consider the free fermion distribution at 
the BCS limit in an asymmetric double-well lattice with 
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FIG. 4: The superlattice configuration to achieve three atoms 
per double well. This is obtained by turning on the lattice 
potentials Vi and V2 simultaneously with relative phase ip > 0, 
producing double wells with a non-zero potential bias between 
the left and right wells. (The overall harmonic potential is 
exaggerated for illustration purposes). In the figure, the solid 
line in each well corresponds to the lowest level, the long 
dotted lines correspond to the Fermi surfaces for "["-atoms and 
|-atoms (with Fermi energies ^ and /ij), which differ due to 
the polarization P > 0, and the dotted rectangles indicate 
those double wells that are occupied by two "[-atoms and one 
|-atom. The /if and /i| are chosen such that j-atoms only 
occupy the left wells while "["-atoms occupy both wells. This 
is the initial configuration needed to measure the hole hopping 
rate t a 3- There is also the possibility of additional "[-atoms 
further from the center of the trap, but the measured molecule 
signal is only sensitive to double wells containing both "[-atoms 
and .[.-atoms. With the conditions given in the text, we insure 
that the only such double wells are those with two "[-atoms 
and one |-atom. 

a bias 612 controlled by the phase shift tp. We would 
like to have two atoms (one spin-| and one spin-|) in 
the deep wells and one spin-f atom in the shallow wells 
as shown in Fig. [U This can be achieved by choosing 
the polarization P and bias 812 so that the atom num- 
bers satisfy iV T > (2\/2 + l)iVo and N l < N , where 

N = (4tt/3) (5 12 /2muj 2 a 2 ) 3/2 . These relations were de- 
rived by requiring that N-\ be great enough that every 
double well which contains a |-atom must also contain 
at least two "[-atoms (so that the molecule signal corre- 
sponds only to double wells containing three atoms), and 
also requiring that no double well contain more than one 
!-atom. (Note that it is not sufficient to require a po- 
larization P > 1/3, since this could be achieved with an 



inner core of double wells containing one "["-atom and one 
!-atom, surrounded by a shell of double wells contain- 
ing only an "[-atom.) jVj must also be small enough that 
there are no more than two j atoms per double well in 
the center of the trap. This condition can be met along 
with the above conditions provided the band gap of the 
lattice is sufficiently great. 

A key assumption in deriving the Hamiltonian ([1]) is 
that in the strongly interacting region there is a signif- 
icant energy gap (of the order of the band gap) which 
separates the four low energy states on each site from 
the other higher energy states Q. With the super- 
lattice technique, one can directly test this assumption 
and measure the energy gap. Given this energy gap, if 
we fill each site with two atoms, there will be no dynam- 
ics as long as the atomic tunneling rate between the two 
sites is small compared with the band gap energy. To 
fill each site with two atoms (one spin-"[, one spin- J,), we 
can start with the free-fermion distribution in the BCS 
limit, choosing the polarization P w and total atom 
number N < AN max (with N max as defined in Eq. 
We then adiabatically turn on V\ and V2 simultaneously 
while keeping a constant ratio V10/V20 > 1- With this 
filling pattern, we should see no dynamics in the strongly 
interacting region, so the atomic distribution over the two 
sites (which will be mapped to the molecular population 
distribution in the BEC limit) will not change with the 
evolution time t. One can also tilt the double-well lattice 
by tuning the bias #12, and measure what is the critical 
612 to turn on the two-site dynamics in the population 
distribution. The measured critical 612 will give an esti- 
mate of the energy gap to excite the system to the high 
energy states. 

In summary, we have described a scheme to test in a 
controllable fashion the predictions of an effective Hamil- 
tonian for strongly interacting fermions in an optical lat- 
tice. With the superlattice technique, one can directly 
test the key assumption in derivation of the Hamiltonian, 
and can measure the physical parameters to confirm the 
particle-assisted tunneling. This scheme provides a quan- 
titative testbed to compare theory with experiments in 
the strongly interacting region. 
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